HerwigH — |- 



Martyn Gigg and Peter Richardson 
Institute of Particle Physics Phenomenology, Department of Physics, University of Durham, Durham, 

DHl 3LE, U.K.. 

We describe the recent development of the Herwig++ event generator. 



Keywords: Hadron Colliders, Monte Carlo Simulations 

1 Introduction 

Monte Carlo event generators have become an essential part of all experimental analyses in 
particle physics. While it has been possible to extend and improve the existing HERWIG^ 
program over many years making further improvements is increasingly difficult. In order to 
include our improved understanding of the physics and recent theoretical developments a pro- 
gramme is therefore underway to produce a new simulation in C++, Herwig++ 12121^ based on 
the same physics philosophy and models. This is part of a wider program within the Monte 
Carlo community to produce a new generation of simulations for the LHC.^ 

We will present the recent physics developments in the Herwig++ simulation, concentrating 
on improvements to the simulation of QCD radiation and new physics, and plans for further 
improvements. 

2 Simulation of QCD Radiation 

The main change between the HERWIG and Herwig++ programs is in the simulation of per- 
turbative QCD radiation. While both programs use an angular-or dere d parton shower designed 
to treat soft-gluon interference effects the algorithm in Herwig++^2EI j^^g ^ number of improve- 
ments: invariance under boosts along the jet axis; improved treatment of radiation from heavy 
quarks; and better coverage of the soft region of phase space. 
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Figure 1; Mass distribution of the quark and lepton in the decay ql —* gX2 ~* 1^^^^ for (a) l'^ and (b) 



In particular the new algorithm uses an improved evolution variable and the quasi-collinear 
splitting functions^ to give better treatment of the radiation from massive particles. In FOR- 
TRAN HERWIG the "dead-cone" approximation of forbidding radiation with angle less than 
m/E, where m is the mass and E the energy of the heavy particle, was used. In Herwig-|--|- this 
is replaced by a smooth suppression of radiation in the direction of the particle. 

The improved treatment of the kinematics of the branchings in the shower means that the 
soft-region of phase space in e~^e~ — > qq is smoothly covered with radiation from the quark and 
anti-quark filling separate regions of phase space which cover the whole region for soft emission 
without overlapping, as was the case with the FORTRAN algorithm. 

A major new feature is the inclusion of radiation from the decaying particle in heavy particle 
decays, for example t — > bW^. This means that in these decays the soft region for gluon emission 
is fully covered, whereas in the FORTRAN program, which only included radiation from the 
decay products, part of the soft region was not filled. This makes correcting the parton shower 
using the exact single emission matrix element simpler. These corrections are now included for 
e+e~ qq, top decay^ and the Drell-Yan process. 

Another key feature of the new algorithm is that momentum reshuffling needed to ensure 
energy and momentum conservation is under greater analytic control which will make it easier to 
match with higher order matrix elements. A number of developments in this area are underway.^ 



3 BSM Physics 

The existing HERWIG program includes a detailed simulation of super symmetric (SUSY) mod- 
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including both spin correlation effects^^ and R-parity violating models.^^^l However, while 
the simulation of SUSY models was highly sophisticated, extending the simulation to other 
models of new physics was difficult and time consuming. In the new simulation we have adopted 
an entirely different approach for the inclusion of new physics models. In the FORTRAN 
simulation the matrix element for each new scattering process and d ecay was added by hand. In 
the new simulation we have included a library based on the HELAS^^ formalism which is used 
in all matrix element calculations. The spin structures for the possible 2^2 matrix elements 
and 1 — *■ 2 decays are included, based on the possible Feynman diagrams for each combina- 
tion of the spins of external particles. The possible scattering processes and decays are then 




Figure 2: Mass distribution of the quark and lepton in the decay g£ qZ* — > l ^l*^ in the UED model for (a) 
and (b) i~ . The results of Herwig++ are compared with the analytic results.'— ^ The mass is given in terms 

of the maximum possible value. 



automatically calculated from the Feynman rules implemented in the code. Using the HELAS 
formalism allows us to include spin correlations in the decays of the fundamental particles, and 
also using the new simulation of tau and hadron decays the correlations in these decays which 
can be important in the decay of SUSY particles and the Higgs boson. 

This approach was originally tested using the Randall- Sundrum model^^ and the Minimal 
Supersymmetric Standard Model (MSSM). An example of the mass distribution of the quark 
and lepton produced in the decay —* qx2 ~^ Q'^^^r; is shown in Fig.[Tl It is important that 
the correlations in this decay are correct as it may be possible to measure the spins of the SUSY 
particles using this decay mode.^^ 

An important test of our new approach for the simulation of BSM physics is the inclusion of 
additional models. We have therefore included the Universal Extra Dimensions model. This 
model is a useful "straw-man" as its particle content is similar to the MSSM but the new particles 
have the same spin as their Standard Model counterparts, rather than opposite spin statistics 
as in SUSY. The mass distribution of the quark and lepton pair in the equivalent decay chain 
to the SUSY chain shown in Fig. [1] is shown in Fig. [2] and compared with the analytic results^^^l 
for these distributions. 



4 Conclusions 



The Herwig++ simulation includes a number of improvements over the previous FORTRAN 
version for both the simulation of perturbative QCD radiation and simulation of new physics. In 
addition to the improvements in the simulation of the perturbative stages of the event generation 
process presented here improvements have been made to the simulation of QED radiation'^ (an 
example of the radiation in Z decays is shown in Fig. [3]) and hadron and tau decays. 

The most recent version of Herwig+-|-'^ is now ready for the simulation of hadron collisions 
and further improvements will be available in the near future. 
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Figure 3: The total energy of the photons radiated in Z boson decays to leptons: (a) shows the spectrum when 
only soft radiation is included and (b) shows the effect of including the coUinear approximation for the hardest 
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